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Organic semiconducting polymers are currently of broad interest as potential low-cost materials 
for photovoltaic and light-emitting display applications. I will give an overview of our work in 
developing a consistent quantum dynamical picture of the excited state dynamics underlying the 
photo-physics. We will also focus upon the quantum relaxation and reogranization dynamics that 
occur upon photoexcitation of a couple of type II donor-acceptor polymer heterojunction systems. 
Our results stress the significance of vibrational relaxation in the state-to- state relaxation and the 
impact of curve crossing between charge- transfer and excitonic states. Furthermore, while a tightly 
bound charge-transfer state (exciplex) remain the lowest excited state, we show that the regeneration 
of the optically active lowest excitonic state in TFB:F8BT is possible via the existence of a steady- 
state involving the bulk charge-transfer state. Finally, we will discuss ramifications of these results 
to recent experimental studied and the fabrication of efficient polymer LED and photovoltaics. 



I. INTRODUCTION 

Organic semiconducting materials have attracted the 
attention of the electronic industry. While the low charge 
mobility compared to crystalline silicon and various other 
characteristics limit the range of potential application for 
organics, a wide range of properties such as solubility in 
organic solvents and the color of light emission can be 
finely tuned using ordinary chemical synthesis. Molecules 
and their arrangement can be easily engineered to fit a 
specific requirement. Moreover, polymer coatings can be 
easily applied to a wide range of surfaces and media, in- 
cluding mechanically flexible ones, to fabricate large scale 
displays or electronic sensors with complex shapes. Or- 
ganic electronics complement and extend the range of 
traditional silicon based electronics. The chemical struc- 
tures of various common semiconducting polymer chains 
are given in Fig^ In practice, PPV is often functional- 
ized with alkyl side chains in order to tune its electronic 
properties and increase its solubility in organic solvents. 

Organic LED devices are typically layered structures 
with luminescent media sandwiched between cathode and 
anode materials which are selected such that their Fermi 
energies roughly match the conduction and valence bands 
of the luminescent material. Often the semiconducting 
media itself consists of a hole transport layer and an 
electron transport layer engineered to facilitate the rapid 
diffusion of the injected carriers away from their image 
charges on the cathode or anode. These carriers are best 
described as polarons since electron-phonon coupling pro- 
duces significant lattice reorganization about the carriers. 
Finally, a third, luminescent layer can be sandwiched be- 
tween the transport layers. In this layer, the electron and 
hole polarons interact and combine to produce excitons. 
The individual spins of the electrons and holes are uncor- 
rected and only singlet excitons are radiatively coupled 
the ground electronic state. In the absence of singlet- 
triplet coupling, this places a theoretical upper-limit or 
1:4 or 25% on the overall efficiency of a LED device and 
it has been long debated whether or not the efficiency of 
organic LED devices is in fact limited by this theoretical 
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FIG. 1: Structures and common short-hand names of various 
conjugated poly-phenylene derived semiconducting polymers 
that are of interest for fabricating luminescent devices. 

upper-limit. 

The electronic properties of these materials are derived 
from the delocalized 7r orbitals found in conjugated poly- 
mers. The 7r electron system is primarily an intramolec- 
ular network extending along the polymer chain. For a 
linear chain, the valance and conduction it and tt* bands 
are typically 1-3 eV wide compared to the intermolccular 
bandwidth (due to 7r-stacking) of about O.leV for well or- 
dered materials. Thus, intrachain charge transport is ex- 
tremely efficient; however, interchain transport typically 
limits the charge mobility for the usual size range of de- 
vices. The polymer backbone is held together through a a 
bonding network. These bonds are considerably stronger 
than the 7r bonds and keep the molecule intact even fol- 
lowing photoexcitation. Hence, we can consider the elec- 
tronic dynamics as taking place within the 7r band and 
treat the localized a bonds as skeletal framework. 

Since the dielectric constant of organic semiconductors 
is relatively low, screening between charges is relatively 
week. At a given radius, r c , thermal fluctuation will be 
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insufficient to break apart an electron/hole pair, 

e 2 

kT= — 

er c 

at 300K, this radius is approximately 20nm, which is on 
the order of a few molecular lengths. If we consider the 
electron/hole pair to be a hydrogenic-type system with 
effective masses equal to the free electron mass and di- 
electric constant of 3, the resulting binding energy is 
about 0.75eV with an effective Bohr radius of 0.3nm, 
which effectively confines the exciton to a single molec- 
ular unit. Finally, if we consider the electron/hole pair 
to be pair of bound Fcrmions, exchange energy result- 
ing from the anti-symetrization of the electron/hole wave 
function splits the spin-singlet and spin-triplet excitons 
by about 0.5-0.7 eV with the spin-triplet lying lower in 
energy than the singlet. While both species are rela- 
tively localized, singlets typically span about lOnm in 
well ordered materials while triplets are much more lo- 
calized. In absence of spin-orbit coupling, emission from 
the triplet states is forbidden. Hence, triplet formation 
in electron/hole capture can dramatically limit the effi- 
ciency of a light-emitting diode device, although strong 
theoretical and experimental evidence indicates that sin- 
glet formation can be enhanced in long-chain polymers. 

Experiments by various groups suggest that in long- 
chain conjugated polymer systems, the singlet exciton 
population can be greatly enhanced and that efficiencies 
as high as 60 to 80% can be easily achieved in PPV type 
systems. [l| On the other hand, in small oligomers, the 
theoretical upper limit appears to hold true. These initial 
experiments were then followed by a remarkable set of ob- 
servations by Wohlgenannt and Vardeny that indicate 
that the singlet to triplet ratio, r > 1 for wide range of 
conjugated polymer systems and that r scales universally 
with the polaron energy-which itself scales inversely with 
the persistence length of the 7r-conjugation, 

r oc 1/n. (1) 

The electro-luminescent efficiency, (f>, is proportional to 
the actual singlet population and is related to r via 

= r/(r + 3). 

Various mechanisms favoring the formation of singlets 
have been proposed for both interchain and intrachain 
e-h collisions. Using Fermi's golden rule, Shuai, Bredas 
et al. indicate that the S cross section for inter- 

chain recombination can be higher than the triplet one 
due to bond-charge correlations. Wohlgenannt et aL[6j 
employ a similar model of two parallel polyene chains. 
Both of these works neglect vibronic and relaxation ef- 
fects. In simulating the intrachain collision of opposite 
polarons, Kobrak and Bittner 0, S @ show that forma- 
tion of singlets are enhanced by the near-resonance with 
the free e-h pair. The result reflects the fact that spin- 
exchange renders the triplet more tightly bound than the 



singlet and hence more electronic energy must be dissi- 
pated by the phonons in the formation of the former. 
The energy-conservation constraints in spin-dependent 
e-h recombination have been analyzed by Burin and 
Ratner|lQ| in an essential-state model. The authors point 
out that nonradiative processes (internal conversion, in- 
tersystem crossing) must entail C=C stretching vibrons 
since these modes couple most strongly to ir — > n* ex- 
citations. Tandon et al. suggest that irrespective of the 
recombination process, interchain or intrachain, the di- 
rect transition to form singlets should always be easier 
than triplets due to its smaller binding energy relative 
to the triplet. A comprehensive review of detailing 
the experiments and theory of this effect was presented 
by Wohlgenannt and Vardenv. [l2l| . By and large, recent 
theoretical models point towards the role of multi-phonon 
relaxation and the scaling of the singlet/triplet splitting 
with chain len gth as dominant factors in determining this 
enhancement. [13, Q HE 

If we assume that the electron/hole capture proceeds 
via a series of microstates one can show that the ratio of 
the singlet to triplet capture cross-sections, r scale with 
the ratio of the exciton binding energies [l7j 

c T 
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If we take the singlet-triplet energy difference to be equal 
to twice the electron- hole exchange energy, AE$t = 2if , 
and expand K in terms of the inverse conjugation length, 

K = K 00 + K^/n+--- 

where 2K oo is the singlet-triplet splitting of an infinitely 
long polymer chain, one obtains 

r oc 1 + nKoo/eg +■■■ 

Since both > 0, r > 1. Moreover, if we take 

eg oc 1/n, we obtain a simple and universal scaling law 
for the singlet-triplet capture ratio, r oc n. This "univer- 
sal" scaling law for r was reported by Wohlgennant and 
Vardeny 0|. What is even more surprising, is that the the 
same scaling law (i.e. slope and intercept for r = an + c) 
describes nearly all organic conjugated polymer systems. 
Hence, it appears that a common set of electronic inter- 
action parameters is transferable between a wide range 
of organic conjugated polymer systems. 

Another general consequence of localized electronic 
states in molecular semiconductors is their effect on the 
molecule itself. Promoting an electron from a ir bonding 
orbital to a ir* anti-bonding orbital decreases the bond 
order over several carbon-carbon bonds. This leads to a 
significant rearrangement of the bond-lengths to accom- 
modate the changes in the electronic structure. By in 
large, for polymers containing phenyl rings, it is the C=C 
bond stretching modes and much lower frequency pheny- 
lene torsional modes that play significant roles in the lat- 
tice reorganization following optical excitation. This is 
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tion threshold as seen by comparing the relative band 
off-sets in FigH Notice that the off-set for TFB:F8BT is 
only slightly larger than 0.5eV, which is approximately 
the exciton binding energy where as in PFB:F8BT the 
off-set is > 0.5eV. Since such blends lie close to the 
stabilization threshold, they are excellent candidates for 
studying the relation between the energetics and the 
kinetics of exciton fission. 

In this paper we will present an overview of our recent 
theoretical work aimed at understanding and modeling 
the state-to-state photophysical pathways in blended het- 
crojunction materials. [lSl Il9| 



FIG. 2: Relative placement of the HOMO and LUMO levels 
for various conjugated polymers. 



evidenced in the strong vibronic features observed in the 
absorption and emission spectra of these materials. 

Finally, one can fabricate devices using blends of semi- 
conducting polymers which phase segregate. For ex- 
ample, the phenylene backbone in F8BT is very pla- 
nar molecule facilitating very delocalized 7r-states. On 
the other hand, TFB and PFB are very globular poly- 
mers due to the tri-amide group in the chain. Conse- 
quently, phase segregation occurs due to more favorable 
7r-stacking interactions between F8BT chains than be- 
tween F8BT and TFB or PFB. Moreover, the electronic 
states in TFB and PFB arc punctuated by the tri-amides. 
This difference in electronic states results in a band off- 
set between the two semiconducting phases. When we 
place the materials in contact with each other, a p-n het- 
erojunction forms. 

One can think of the HOMO and LUMO energy lev- 
els of a given polymer as corresponding to the top and 
bottom of the valance and conduction bands respectively. 
For the polymers under consideration herein, the relative 
band edges are shown in Fig. [21 In Type II heterojunc- 
tion materials, the energy bands of the two materials are 
off set by AE. If the exciton binding energy eb > AE, 
excitonic states will the lowest lying excited state species, 
resulting in a luminescent material with the majority of 
the photons originating from the side with the lowest op- 
tical gap. Since the majority of the charge carriers are 
consumed by photon production, very little photocurrent 
will be observed. 

On the other hand, charge transfer states across the 
interface will be energetically favored if Eb < AE. Here, 
any exciton formed will rapidly decay in to a charge- 
separated state with the electron and hole localizing on 
either side of the junction. This will result in very little 
luminescence but high photocurrent. Consequently, het- 
erojunctions of PPV and BBL which have a large band 
off-set relative to the exciton binding energy are excellent 
candidate materials for organic polymer solar cells. 

Heterojunctions composed of TFB:F8BT and 
PFB:F8BT lie much closer to the exciton stabiliza- 



II. TWO-BAND CONFIGURATION 
INTERACTION MODEL 

Our basic description is derived starting from a model 
for the on-chain electronic excitations of a single conju- 
gated polymer chain. 0,H3,|^| This model accounts for 
the coupling of excitations within the 7r-orbitals of a con- 
jugated polymer to the lattice phonons using localized 
valence and conduction band Wannicr functions (\h) and 
\p)) to describe the it orbitals and two optical phonon 
branches to describe the bond stretches and torsions of 
the the polymer skeleton. 

H = E^-n + V mn )AlA a 

mn 

nmi/( x ^ 
+ J^^l^ + Kli^i+l^) +Pin ( 2 ) 

where A\ and A n are Fermion operators that act upon 
the ground electronic state |0) to create and destroy elec- 
tron/hole configurations \n) — \hp) with positive hole in 
the valence band Wannicr function localized at h and an 
electron in the conduction band Wannier function p. Fi- 
nally, qi^ and p^ correspond to lattice distortions and 
momentum components in the «-th site and ^-th optical 
phonon branch. 

Wannier functions are essentially spatially localized 
basis functions that can be derived from the band- 
structure of an extended system. Quantities such as 
the exchange interaction and Coulomb interaction can be 
easily computed within the atomic orbital basis; however, 
there are many known difficulties in computing these 
within the crystal momentum representation. Because 
of this, is is desirable to develop a set of orthonormal 
spatially localized functions that can be characterized by 
a band index and a lattice site vector, R^. These are the 
Wannier functions, which we shall denote by a n (r — R^) 
and define in terms of the Bloch functions 

a ^ r - R ^ = tJ^TI / e- lkR ^ nk (r)dk. (3) 
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The integral is over the Brillouin zone with volume 
V = (2n) d /^l and fl is the volume of the unit cell (with 
d dimensions). A given Wannier function is defined for 
each band and for each unit cell. If the unit cell happens 
to contain multiple atoms, the Wannier function may be 
delocalized over multiple atoms. The functions are or- 
thogonal and complete. 



The Wannier functions are not energy eigenfunctions of 
the Hamiltonian. They are, however, linear combinations 
of the Bloch functions with different wave vectors and 
therefore different energies. For a perfect crystal, the 
matrix elements of H in terms of the Wannier functions 
are given by 



J 'af(r - R v )H a n (r - RJdr = J e^ R -- kR ^^ lk (r)H ^ nk (r)drdqdk 

= £ n {R v - R„)6 nl (4) 

where 

£ n {R v -R„) = J e ik ^- R ^E n (k)dk. 

Consequently, the Hamiltonian matrix elements in the Wannier representation are related to the Fourier components 
of the band structure, E n (k). Therefore, given a band structure, we can derive the Wannier functions and the single 
particle matrix elements, F^ n . 



The single-particle terms, i^ n , are derived at the 
ground-state equilibrium configuration, = 0, from the 
Fourier components f r and f r of the band energies in 
pseudomomentum space. 

F mn = $mti(m\f\n) - 8 mn (m\f\n) 

^mnfrn—n ^Tnnf^n—Ji (5) 

Here, f mn and f mn are the localized energy levels 
and transfer integrals for conduction-band electrons and 
valence-band holes. At the ground-state equilibrium ge- 
ometry, = 0, these terms can be computed as Fourier 
components of the one-particle energies in the Brillouin 
zone. For example, for the conduction band: 

fmn = fm-n = ^ [ e k ^ m ~^ dk (6) 
o z J Bz 

where k is the pseudomomentum for a 1-dimcnsional 
lattice with unit period. For the case of cosine-shaped 
bands, 

e(fc) = f + 2A cos(fc) 

the site-energies are given by the center f a and the trans- 
fer integral between adjacent Wannier functions is given 
by fx. The band-structure and corresponding Wannier 
functions for the valence and conduction bands for PPV 
are shown in Fig. El ii El El For the intrachain terms, 
we use the hopping terms and site energies derived for 
isolated polymer chains of a given species, t^u, where 
our notation denotes the parallel hopping term for the 
ith chain (i = 1,2). For PPV and similar conjugated 
polymer species, these are approximately 0.5eV for both 
valence and conduction 7r bands. 
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FIG. 3: Computed band-structure and vinylene-centered 
Wannier functions for PPV. 



The two-particle interactions are spin dependent with 



VZi = -( mn \\ nm ) (7) 
V„ S ln = V^ n + 2(mn\\mn). (8) 
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TABLE I: Electron/hole integrals for poly-phenylene-type 
polymer chains. 



Term 


Functional form 


Parameters 


Direct Coloumb 


J(r) = Jo/(l + r/r a ) 


J = 3.092eV 






r = 0.6840a 


Exchange 


K(r) = K e~ r/r ° 


K = 1.0573eV 






r a = 0.4743a 


Dipole-Dipole 


D(r) = D (r/r )- 3 


D = -0.03209eV 






r = 1.0a 



Note: a = unit lattice spacing. 



for triplet and singlet combinations respectively with 

(mn\\ij) = JdlJ d2C(l)^(2) u (12)0 l (l)0-(2) 

With the exception of geminate WFs, orbital overlap 
is small such that the two-body interactions are lim- 
ited to Coulomb, J(r) and exchange integrals, K(r) re- 
flecting e-h attraction and spin-exchange coupling non- 
geminate configurations, and transition dipolc-dipole in- 
tegrals D(r) coupling only geminate singlet electron-hole 
pairs. TableQlgives a listing of the electron-hole integrals 
and their parameters we have determined for PPV and 
similar poly-phenylene based conjugated chains. We have 
found that these are quite transferable amongst this class 
of conjugated polymers and allow us to focus upon mod- 
eling similar poly-phenylene chains through variation of 
the Wannier function band-centers (i.e. site-energies) 
and band- widths (i.e. intrachain hopping integrals). 

Since we will be dealing with inter-chain couplings, 
we make the following set of assumptions. First, the 
single-particle coupling between chains is expected to 
be small compared to the intramolecular coupling. For 
this, we assume that the perpendicular hopping inte- 
gral t± = 0.01 eV. This is consistent with LDA calcu- 
lations performed by Vogl and Campbell and with the 
t_i_ ~ 0.15/i estimate used in an earlier study of inter- 
chain excitons by Yu et al. |22L |23[ Furthermore, we as- 
sume that the J(r), K(r), and D(r) two-particle interac- 
tions depend only upon the linear distance between two 
sites, as in the intrachain case. Since these are expected 
to be weak given that the interchain separation, d, is 
taken to be some what greater than the inter-monomer 
separation. 

An important component in our model is the coupling 
between the electronic and nuclear degrees of freedom. 
These we introduce via a linear coupling term: 



dfn 



f(2^ 3 ) 1/2 (* 



(9) 



where S is the Huang- Rhys factor which can be obtained 
from vibronic features in the experimental photoemis- 
sion spectra. For the case of conjugated polymers such 



as PPV and similar poly-phenylene species, the emis- 
sion spectra largely consists of a series of well-resolved 
vibronic features corresponding to the C=C stretching 
modes in the phenylene rings with typical Huang Rhys 
factors of 0.4. In addition, low frequency vibrational 
modes (due to torsions) contribute a broad featureless 
background. Consequently, we include two intramolecu- 
lar optical phonon branches which correspond roughly to 
the high-frequency C=C bond stretching modes within 
a given repeat unit and a second low- frequency mode, 
which in the case of PPV are taken to represent the 
phenylene torsional modes. The electron-phonon cou- 
plings are assumed to be transferable between the vari- 
ous chemical species. Since the modes are assumed to be 
intramolecular, we do not include interchain couplings in 
the phonon Hessian matrix. 

Upon transforming H into the diabatic representation 
by diagonalizing the electronic terms at qi^ — 0, we ob- 
tain a series of vertical excited states \a a ) with energies, 
e° and normal modes, with frequencies, lj^. (We will 
assume that the sum over £ spans all phonon branches). 



H = J2e° a \a )(a \ + g° ab( qd\ a o) (bo\ + \b )(a \) 

a ab£ 

+ JEK 2 ^ + P D- ( 10 ) 



The adiabatic or relaxed states can be determined then 
by iteratively minimizing e a (Q^) — (a\H\a) according to 
the self-consistent equations 



9aai + ^iQi = 0. 



(11) 



Thus, each diabatic potential surface for the nuclear 
lattice motion is given by 



(12) 



These are shown schematically as S a and Sb in Fig.|3|with 
being a collective normal mode coordinate. On can 
also view this figure as a slice through an N— dimensional 
coordinate space along normal coordinate In this fig- 
ure, e° and eg are the vertical energies taken at the 
ground-state equilibrium geometry = 0. The adia- 
batic energies, taken at the equilibrium geometry of each 
excited state are denoted as e a and eb- While our model 
accounts for the distortions in the lattice due to elec- 
tron/phonon coupling, we do not account for any adi- 
abatic change in the phonon force constants within the 
excited states. Lastly, the electronic coupling between 
diabatic curves is given by g ab which we compute at 
the ground-state geometry (g° b ). We assume that both 
the diagonal g° a and off-diagonal g° b terms can be de- 
rived from the spectroscopic Huang-Rhys parameters as 
in Eq. H 

The advantages of our approach is that it allows us to 
easily consider the singly excited states of relatively large 
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FIG. 4: Schematic representation of excited state Diabatic 
potentials obtained within our approach. The ground state 
configuration is taken as — with vertical excitation ener- 
gies ate£ and e£ and adiabatic (minimum) energies at e a and 



conjugated polymer systems. Our model is built from 
both ah initio and experimental considerations and can 
in fact reproduce most of the salient features of the vi- 
bronic absorption and emission spectra for these systems. 
The model is limited in that we cannot include specific 
chemical configurational information about the polymers 
other than their conjugation length and gross topology. 
For isolated single chains, the model is rigorous. For 
multiple chains, our interchain parameterization does not 
stand on such firm ground since technically the Wan- 
nier functions are derived from a quasi-one-dimensional 
band structure. Nonetheless, our model and results pro- 
vide a starting point for predicting and interpreting the 
complex photophysical processes within these systems. 
We next move on to describing the state-to-state inter- 
conversion proceses that occur following both photo and 
electro-excitation. 



jugated polymers in which both electroluminescence and 
photolumincsccncc originates from the same S\ — S x t 
state. This implies that internal conversion dynamics 
are fast relative to the fluorescence lifetime. 

Coupling the electronic relaxation dynamics to the vi- 
brational dynamics is a formidable task. An exact quan- 
tum mechanical description of this is currently well be- 
yond the state of the art of current computational meth- 
ods. One can, however, compute the state-to-state rate 
constants using Fermi's golden rule and arrive at a rea- 
sonable picture. 

If we assume that the vibrational bath described by 
H p h remains at its ground-state geometry, then the state- 
to-state transition rates are easily given by Fermi's golden 
rule: 



7r£ff(l 



n(w a6 ))(r(w 4 - uj ab ) - T(w£ + w a b)\13) 



where n(u>) is the Bose-Einstein population for the 
phonons, T is a Lorentzian broadening in which the 
width is inversely proportional to the phonon lifetime 
used to smooth the otherwise discrete phonon spectrum, 
Lu a i, = (e° — e^)/h. In order for a transition to occur, there 
must be a phonon of commensurate energy to accommo- 
date the energy transfer. The coupling term, g a b£, is the 
diabatic coupling in the diabatic Hamiltonian given in 

Eq.uni 

This static model is fine so long as either the nuclear 
relaxation has little effect on the state to state rate con- 
stant or if the electronic transitions occur on a time-scale 
which is short compared to the nuclear motion. However, 
if lattice reorganization does play a significant role, then 
we need to consider the explicit nuclear dynamics when 
computing the state-to-state rates. If we assume that vi- 
brational relaxation within a given diabatic state is rapid 
compared to the inter-state transition rate, we can con- 
sider the transitions as occurring between displaced har- 
monic wells 



k a b = ^j- \V ab \ 2 F 



(14) 



where V a b is the coupling between electronic states a and 
b and 



III. STATE TO STATE RELAXATION 
DYNAMICS 



x S(e a (v a ) - e{v b ) + AE ab ) (15) 



The electronic levels in our model are coupled to the 
lattice phonons as well as the radiation field. Conse- 
quently, relaxation from a given electronic state can oc- 
cur via state to state interconversion via phonon excita- 
tion or absorption or fluorescent decay to the S ground 
state. For triplet excitations, only phonon transitions are 
allowed. For the singlets, fluorescence occurs primarily 
from the lowest S n state independent of how the exci- 
tation was prepared. This certainly holds true for con- 



is the thermally averaged Franck-Condon weighted den- 
sity of nuclear vibrational states. Here, v a and Vb denote 
the vibronic states, Pth is the Boltzmann distribution 
over the initial states, e a (v a ) and e(i>b) are the corre- 
sponding enregies, and AE a b is the electronic energy gap 
between a and b. In the classical limit, T becomes 



HEab) = 



1 



^4TrE x k B T 



exp 



(Ex + AE ab 
4E x k B T 



(16) 
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where E\ is the reorganization energy as sketched in 
Fig.H 

Each of these terms can be easily computed from the 
diabatic Hamiltonian in Eq. I1UI The diabatic coupling 
matrix element between the adiabatically relaxed excited 
states, |V^b| 2 , requires some care since we are considering 
transitions between eigenstates of different Hamitonians 
(corresponding to different nuclear geometries). Since 
the vertical — states provide a common basis, |a°), 
we can write 

V ab = J2(*\*°)9a b (b°\b) (17) 

a°b° 

where g° b is the diabatic matrix element computed at the 
equilibrium geometry of the ground-state. 

Once we have the rate constants computed, it is a sim- 
ple matter to integrate the Pauli master equation for the 
state populations. 

P a (t) = Y / ( k baPb-k ab P a )-k r a ad P a (18) 
b 

where k r a ad is the radiative decay rate of state a. 

tr = g! (1+nKo)) ^ (19) 

where /i o are the transition dipoles of the excited sin- 
glets. These we can compute directly from the Wannier 
functions or empirically from the photoluminescence de- 
cay rates for a given system. Photon mediated transi- 
tions between excited states are highly unlikely due to 
the Lo^ b density factor of the optical field. In essence, so 
long as the non-equilibrium vibrational dynamics is not 
a decisive factor, we can use these equations to trace the 
relaxation of an electronic photo- or charge-transfer ex- 
citation from its creation to its decay including photon 
outflow measured as luminescence. 



IV. EXCITON REGENERATION DYNAMICS 

Donor-acceptor heterojunctions composed of blends of 
TFB with F8BT and PFB with F8BT phase segregate to 
form domains of more or less pure donor and pure accep- 
tor. Even though the polymers appear to be chemically 
quite similar, the presence of the tri-phenyl amine groups 
in TFB and PFB cause the polymer chain to be folded up 
much like a carpenter's rule. F8BT, on the other hand, 
is very rod-like with a radius of gyration being more or 
less equivalent to the length of a give oligomer. Molecu- 
lar dynamics simulations of these materials by our group 
indicate that segregation occurs because of this differ- 
ence in morphology and that the interface between the 
domains is characterized by regions of locally ordered n- 
stacking when F8BT rod-like chains come into contact 
with more globular PFB or TFB chains. 

As discussed earlier, TFB:F8BT and PFB:F8BT sit 
on either side of the exciton destabalization threshold. 
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FIG. 5: Vertical and relaxed energies of the lowest lying states 
in the TFB:F8BT and PFB:F8BT heterojunctions. In each, 
A and C refer to the interchain exciplex state and B and D 
refer to the predominantly intrachain F8BT exciton state. 



In TFB:F8BT, the band off-set is less than the exciton 
binding energy and these materials exhibit excellent LED 
performance. On the other hand, devices fabricated from 
PFB:F8BT where the exciton binding energy is less than 
the off-set, are very poor LEDs but hold considerable 
promise for photovoltic devices. In both of these systems, 
the lowest energy state is assumed to be an interchain 
exciplex as evidenced by a red-shifted emission about 50- 
80ns after the initial photoexcitation. [24| In the case 
of TFB:F8BT, the shift is reported to be 140±20 meV 
and in PFB:F8BT the shift is 360±30 meV relative to 
the exciton emission, which originates from the F8BT 
phase. Bearing this in mind, we systematically varied 
the separation distance between the cofacial chains from 
r = 2a— 5a (where a = unit lattice constant) and set tj_ = 
0.01 in order to tune the Coulomb and exchange coupling 
between the chains and calibrate our parameterization. 

For large interchain separations, the exciton remains 
localized on the F8BT chain in both cases. As the chains 
come into contact, dipole-dipole and direct Coulomb cou- 
plings become significant and we begin to see the effect 
of exciton destabilization. For TFB:F8BT, we select and 
interchain separation of r = 2.8a giving a 104meV split- 
ting between the vertical exciton and the vertical exciplex 
and 87.4 meV for the adiabatic states. For PFB:F8BT, 
we chose r — 3a giving a vertical exciton-exciplex gap 
of 310 meV and an adiabatic gap of 233 meV. In both 
TFB:F8BT and PFB:F8BT, the separation produce in- 
terchain exciplex states as the lowest excitations, with 
energies reasonably close to the experimental shifts. 

Fig . If V I compares the vertical and adiabatic energy lev- 
els in the TFB:F8BT and PFB:F8BT chains and Figs.© 
and show the vertical and relaxed exciton and charge- 
separated states for the two systems. Here, sites 1-10 
correspond to the TFB or PFB chains and 11-20 cor- 
respond to the F8BT chain. The energy levels labeled 
in Fig. IIVI correspond to the states plotted in Figs. 
and0 We shall refer to states A and B as the vertical 
exciplex and vertical exciton and to states C and D as 
the adiabatic exciplex and adiabatic exciton respectively. 
Roughly, speaking a pure exciplex state will have the 
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FIG. 6: Excited state electron/hole densities for TFB:F8BT 
heterojunction. The electron/hole coordinate axes are such 
that sites 1-10 correspond to TFB sites and 11-20 correspond 
to F8BT sites. Note the weak mixing between the interchain 
charge-separated states and the F8BT exciton in each of these 
plots. 



charges completely separated between the chains and will 
contain no geminate electron/hole configurations. Like- 
wise, strictly speaking, a pure excitonic state will be lo- 
calized to a single chain and have only geminate elec- 
tron/hole configurations. Since site energies for the the 
F8BT chain are modulated to reflect to internal charge- 
separation in the F8BT co-polymer as discussed above, 
we take our "exciton" to be the lowest energy state that is 
localized predominantly along the diagonal in the F8BT 
"quadrant" . 

In the TFB:F8BT junction, the lowest excited state 
is the exciplex for both the vertical and adiabatic lat- 
tice configurations with the hole on the TFB and the 
electron on the F8BT. (FigHJY,C) In the vertical case, 
there appears to be very little coupling between intra- 
chain and interchain configurations. However, in the adi- 
abatic cases there is considerable mixing between intra- 
and inter-chain configurations. First, this gives the adia- 
batic exciplex an increased transition dipolc moment to 
the ground state. Secondly, the fact that the adiabatic 
exciton and exciplex states are only 87 meV apart means 
that at 300K, about 4% of the total excited state popu- 
lation will be in the adiabatic exciton. 

For the PFB:F8BT heterojunction, the band off-set is 
greater than the exciton binding energy and sits squarely 
on the other side of the stabilization threshold. Here the 
lowest energy excited state (Figs.[7J\ and B) is the inter- 



FIG. 7: Excited state electron/hole densities for PFB:F8BT 
heterojunction. The axes are as in previous figures except 
that sites 1-10 correspond to PFB sites and sites 11-20 to 
F8BT sites. 



chain charge-separated state with the electron residing on 
the F8BT (sites 11-20 in the density plots in Fig0) and 
the hole on the PFB (sites 1-10). The lowest energy ex- 
citon is almost identical to the exciton in the TFB:F8BT 
case. Remarkably, the relaxed exciton (Fig. [T)D) shows 
slightly more interchain charge-transfer character than 
the vertical exciton (Fig. [7p). While the system read- 
ily absorbs at 2.3eV creating a localized exciton on the 
F8BT, luminescence is entirely quenched since all popu- 
lation within the excited states is readily transfer to the 
lower-lying interchain charge-separated states with van- 
ishing transition moments to the ground state. 

In calculating the state-to-state interconversion rates 
for TFB:F8BT, we note two major differences between 
the static and the adiabatic Marcus-Hush approaches 
(See Fig. [SJ ■ First is the sparsity of the latter with tran- 
sitions being limited to states with smaller energy differ- 
ences. This leads to a relaxation dynamics that is more 
intertwined with the DOS. Second is the relatively faster 
rates calculated in the latter leading to interconversion 
lifetimes in the femtosecond (fs) to a couple of picosecond 
(ps) regime as opposed to hundreds of ps in the former. 
The same general difference is observed for PFB:F8BT 
(Not shown). These marked differences in the distribu- 
tion of rates and their range of magnitudes are brought 
about by the introduction of the reorganization energy 
as a parameter in the rates calculation to complement 
the energy differences between the states. It provides a 
way to incorporate lattice distortions in the semiclassi- 
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FIG. 8: (Color online) TFB-F8BT internal conversion rates 
distribution at 290 K. Rates are in ps -1 . Note the sparsity 

and relatively faster Marcus rates compared to the diabatic FIG. 9: Time-evolved populations of the lowest CT (solid 
rates. lines) and XT (dashed lines) relaxed states of TFB:F8BT(TF) 

and PFB:F8BT(PF) in the Marcus-Hush approach at 290 K. 



cal limit into the relaxation dynamics. While this is not 
fully dynamical in its account of the lattice distortions, it 
improves upon the static approximation previously em- 
ployed. 

The photoexcitation of heterojunction systems is simu- 
lated by populating a higher-lying excitonic state. Fig. 
shows the time-evolved populations of the lowest charge- 
transfer (CT) and excitonic (XT) vertical and relaxed 
states, respectively, in photoexcited TFB:F8BT and 
PFB:F8BT. We see, in both approaches, that the relax- 
ation to the lowest CT state is faster in TFB:F8BT than 
in PFB:F8BT. Furthermore, the relaxation from the XT 
state to the CT state occurs faster in the former. This is 
despite the XT state being formed faster in the latter for 
both cases. This is manifested more in the Marcus-Hush 
approach shown in Fig. where despite reaching a maxi- 
mum population of 0.86 in 250 fs as opposed to just 0.40 
in 500 fs, the XT— >CT interconversion is practically done 
in 2 ps in TFB:F8BT compared to 10 ps in PFB:F8BT. 
In addition, we note that the XT state reaches a steady 
state population in TFB:F8BT whereas it goes to zero 
in PFB:F8BT. This small but non-zero population of the 
XT state is consistent with the distributed thermal pop- 
ulation of states of 0.022 at 290 K owing to the fact that 
this XT state is 95 meV higher in energy relative to the 
lowest CT state, pi 

Interestingly, while the overall XT— >CT interconver- 
sion occurs in just a couple of ps in both heterojunction 
systems, a closer look into the rates reveal that this relax- 
ation does not occur directly. Rather, it involves the next 
lowest CT state. Fig. ^] show the relevant interconver- 
sions between the three lowest states of both systems: the 
lowest CT state(CTl), the next lowest CT state(CT2), 
and the lowest XT state (XT). It is worth noting that the 



considerable mixing between the intra-chain and inter- 
chain configurations of the former compared to those of 
the latter. In TFB:F8BT, the direct XT^CTl transition 
(~ 10~ 3 ps -1 ) is at least 3 orders of magnitude slower 
than the corresponding XT^CT2^CT1 transition route 
(> 1 ps -1 ), the indirect route being consistent with the 
evolution data (Fig. EJ. Thus, the XT— >CT conversion 
occurs via the CT2 state and not directly. The reverse 
transitions for both routes are slower but have the same 
order of magnitude difference. This CT1^CT2^XT 
transition (~ 10 _1 ps -1 ) effectively presents a regener- 
ation pathway for the XT. This leads to an XT state 
population that is always at equilibrium with the CT1 
state. 

In PFB:F8BT the XT^CTl and XT^CT2 conversion 
occur at relatively the same rate (~ 10 _1 ps -1 ) while 
their reverse transitions are at least 2 orders of magnitude 
slower. Consequently, XT is not regenerated. The role 
played by CT2 as a bridge state is apparently relative 
to whether it has a slightly higher or lower energy than 
the XT as has been accounted by Morteani et aZ. |24ll25j . 
Spontaneous transition rates are typically faster when 
going from a higher to a lower energy state than the 
reverse according to detailed balance. Here, CT1 is the 
exciplex state which exhibit sizable mixing with the bulk 
CT state (CT2). When CT2 has a higher energy than 
XT, such as in TFB:F8BT, a fraction of the population 
in CT2 converts to XT. If it has a lower energy relative 
to the XT state such as in PFB:F8BT, this regeneration 
of the XT, practically, does not occur. 

To see the effect of temperature, the interconversion 
rates were calculated at 230, 290, and 340 K. Fig. El 
shows how the interconversions among the three lowest 
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FIG. 10: Relevant Marcus-Hush interconversion rates for the 
3 lowest states of (left) TFB:F8BT and (right) PFB:F8BT. 
In both cases, relaxation proceeds from the density of states 
(DOS) to the lowest excitonic state(XT) (off-set to the right 
relative to the CT states for clarity of relaxation route) be- 
fore relaxing to the lowest charge-transfer state(CTl). CT1 
proceeds to equilibrium with the next higher CT state(CT2). 
In PFB:F8BT, CT2 has a lower energy than XT where as 
in TFB:F8BT, it has a higher energy. Also shown are the 
radiative rates emanating from the XT states which are 
strongly coupled to ground state, So, of both systems and 
the TFB:F8BT CT1 state which is just weakly coupled to So. 



states (CT1, CT2, and XT) of TFB:F8BT, as illustrated 
in Fig. 1101 vary with temperature. This dependence is 
given in an Arrhenius plot of log k vs 1/T and gives a 
linear plot for each transition having a slope associated 
with the activation energy, E ac t , for that particular tran- 
sition. This activation energy has the expression, 



{AE- 



4£a 



(20) 



Transitions to lower energy states are given as solid lines 
while those going to higher energy states are given as 
dashed lines. Curiously, although XT^CTl is exother- 
mic compared to XT^CT2 which is endothermic, the 
latter is a more favorable transition. This has to do with 
the fact that XT^CTl has an activation energy almost 
three times greater than that of XT— >CT2. As alluded 
to above, this is a consequence of the former being in the 
inverted region while the latter being in the normal re- 
gion. In the inverted region, the larger AE is, the larger 
E act as opposed to the more familiar normal region where 
E act decreases as AE increases. Having stated this, how- 
ever, we note that in the former, due to maximal overlap 
between the vibrational modes of the two states, transi- 
tions may be possible via tunneling processes. Overall, 
in TFB:F8BT, we see an increase in the fraction of the 
total excited state population in XT as temperature in- 
creases. At 230 K only 0.81% is in XT while at 290 and 
340 K, 2.16% and 3.67% is in XT, respectively. 

Finally, we note that all transition rates increase with 
temperature except for the CT2^XT in TFB:F8BT and 
XT^CT2 in PFB:F8BT which decrease with tempera- 
ture. Such a trend, while not uncommon in chemical 
reactions, are though to be indicative of a more compli- 
cated transition mechanism as noted by Porter [2^. We 
surmise this to be due to the coupling between the low 
frequency vibrational modes of the initial state with the 
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FIG. 11: Interconversion rates between the 3 lowest states 
of TFB:F8BT as a function of temperature (230, 290, and 
340 K). Plot is given as log k vs 1/T. Transitions to lower en- 
ergy states are given as solid lines while their reverse are given 
as dashed lines. The CT2«->CT1, CT2^XT and XT^CTl 
are plotted as squares, triangles and circles, respectively. All 
transition rates increase directly with temperature except 
the CT2^XT conversion which decreases as temperature in- 
creases. 



high frequency vibrational modes of the final state as in 
the case of an early transition state in reactive scatter- 
ing. What we are definite of is that both transitions are 
exothermic and in the normal region which makes one 
speculate if there is any connection. 



V. DISCUSSION 

In this paper, we gave an overview of our recent work 
in developing a theoretical understanding interfacial exci- 
tonic dynamics in a complex material system. The results 
herein corroborate well with the experimental results on 
these systems. In particular, following either charge in- 
jection or photo-excitation, the system rapidly relaxes to 
form the interchain charge-separated species. In the ex- 
perimental data, this occurs within the first 10 or so ps 
for the bulk material. Our calculations of a single pair 
of cofacial chains puts the exciplex formation at about 
1 ps. Likewise, the experimental time-resolved emission 
indicates the regeneration occurs on a much longer time- 
scale with most of the time-integrated emission coming 
from regenerated excitons. This too, is shown in our 
calculations as evidenced by the slow thermal repopu- 
lation of the XT state in the TFB:F8BT system. Since 
this state has a significant transition dipole to the ground 
state, population transfered to this state can either decay 
back to the CT state via thermal fluctuations or decay 
to the ground state via the emission of a photon. Since 
this secondary emission is dependent upon the thermal 
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population of the XT state at any given time, the effi- 
ciency of this process shows a strong dependency upon 
the temperature of the system. 

The exciton model we present herein certainly lacks 
the molecular level of details so desired by materials 
chemists. However, it offers a tractable way of build- 
ing from molecular considerations the salient physical in- 
teractions that give rise to the dynamics in the excited 
states of these extended systems. In building this model 
we make a number of key assumptions. First, and per- 
haps foremost, that the excited states are well described 
via bands of 7r orbitals and that from these bands we 
can construct localized Wannier functions. Hand in hand 
with this assumption is that within the general class of 
oligo-phenylene derived polymers, configuration interac- 
tion matrix elements, hopping integrals, elect ron/phonon 
couplings, and phonon spectra are transferable from one 
system to another. This is a fairly dangerous approxi- 
mation since it discounts important contributions from 
heteroatoms, side-chains, and chain morphology. How- 
ever, given that a single oligomeric chain of F8BT with 
10 repeat units has well over 300 atoms, such potentially 
dire approximations are necessary in order to extract the 
important features of these very extended systems. 

Secondly, we make the assumption that the explicit 
vibrational dynamics can be integrated out of the equa- 
tions of motion for the electronic states. This is prob- 
ably not too extreme of an assumption so long as we 



can assume that the phonons remain thcrmalized over 
the course of the electronic relaxation. However, looking 
back at the level correlation diagrams, crossings between 
diabatic states are present in this system and hence con- 
ical intersections between electronic states may play an 
important role. Finally, we discount the effects of elec- 
tronic coherence. This, too, may have a profound im- 
pact upon the final state-to-state rate constants since it 
is well recognized that even a small amount of quantum 
coherence between states leads to a dramatic increase 
in the transition rate. Fortunately, many of the papers 
presented in this proceedings address these assumptions. 
Approaches, such as the MCTDH method presented by 
Thoss, the DFT based non-adiabatic molecular dynam- 
ics approach (NAMD) developed by Prezhdo, for exam- 
ple, are important strides towards achieving a molecular 
level understanding of complex photophysical processes 
in light-emitting and light-harvesting materials. 
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